Detecting pore size distribution of activated carbon by low-field nuclear magnetic resonance

In this study, the transverse relaxation time (T₂) of activated carbon (AC) in different relative environment humidity was detected firstly by low-field nuclear magnetic resonance (LFNMR). The pore size (diameter) of AC distributions was calculated by the relationship between T₂ and surface relaxation rate (ρ), where ρ was obtained by the detection of nine porous materials with known pore size. The results showed that the pore size distributions of AC calculated by ρ < 0.19 nm/ms were in good agreement with that obtained by nitrogen adsorption method and proved that LFNMR as a new detection method was feasible for characterizing AC pore size distribution.     

Effect of wall slip on the viscoelastic particle ordering in a microfluidic channel

The formation of a line of equally spaced particles at the centerline of a microchannel, referred as “particle ordering,” is desired in several microfluidic applications. Recent experiments and simulations highlighted the capability of viscoelastic fluids to form a row of particles characterized by a preferential spacing. When dealing with non-Newtonian fluids in microfluidics, the adherence condition of the liquid at the channel wall may be violated and the liquid can slip over the surface, possibly affecting the ordering efficiency. In this work, we investigate the effect of wall slip on the ordering of particles suspended in a viscoelastic liquid by numerical simulations. The dynamics of a triplet of particles in an infinite cylindrical channel is first addressed by solving the fluid and particle governing equations. The relative velocities computed for the three-particle system are used to predict the dynamics of a train of particles flowing in a long microchannel. The distributions of the interparticle spacing evaluated at different slip coefficients, linear particle concentrations, and distances from the channel inlet show that wall slip slows down the self-assembly mechanism. For strong slipping surfaces, no significant change of the initial microstructure is observed at low particle concentrations, whereas strings of particles in contact form at higher concentrations. The detrimental effect of wall slip on viscoelastic ordering suggests care when designing microdevices, especially in case of hydrophobic surfaces that may enhance the slipping phenomenon.     

Imposing accurate wall boundary conditions in corrective‐matrix‐based moving particle semi‐implicit method for free surface flow

Corrective matrix that is derived to restore consistency of discretization schemes can significantly enhance accuracy for the inside particles in the Moving Particle Semi‐implicit method. In this situation, the error due to free surface and wall boundaries becomes dominant. Based on the recent study on Neumann boundary condition (Matsunaga et al, CMAME, 2020), the corrective matrix schemes in MPS are generalized to straightforwardly and accurately impose Neumann boundary condition. However, the new schemes can still easily trigger instability at free surface because of the biased error caused by the incomplete/biased neighbor support. Therefore, the existing stable schemes based on virtual particles and conservative gradient models are applied to free surface and nearby particles to produce a stable transitional layer at free surface. The new corrective matrix schemes are only applied to the particles under the stable transitional layer for improving the wall boundary conditions. Three numerical examples of free surface flows demonstrate that the proposed method can help to reduce the pressure/velocity fluctuations and hence enhance accuracy further.     

A modified SPH method to model the coalescence of colliding non-Newtonian liquid droplets

This paper develops a modified smoothed particle hydrodynamics (SPH) method to model the coalescence of colliding non-Newtonian liquid droplets. In the present SPH, a van der Waals (vdW) equation of state is particularly used to represent the gas-to-liquid phase transition similar to that of a real fluid. To remove the unphysical behavior of the particle clustering, also known as tensile instability, an optimized particle shifting technique is implemented in the simulations. To validate the numerical method, the formation of a Newtonian vdW droplet is first tested, and it clearly demonstrates that the tensile instability can be effectively removed. The method is then extended to simulate the head-on binary collision of vdW liquid droplets. Both Newtonian and non-Newtonian fluid flows are considered. The effect of Reynolds number on the coalescence process of droplets is analyzed. It is observed that the time up to the completion of the first oscillation period does not always increase as the Reynolds number increases. Results for the off-center binary collision of non-Newtonian vdW liquid droplets are lastly presented. All the results enrich the simulations of the droplet dynamics and deepen understandings of flow physics. Also, the present SPH is able to model the coalescence of colliding non-Newtonian liquid droplets without tensile instability.     

An Insight into the Diffusion of Unsaturated Polyester (UP) Resin into Expanded Graphite (EG) to Improve the Mechanical Properties of the UP/EG Composites

Abstract

Polymer/expanded graphite (EG) nanocomposites have great importance in many industrial applications mainly due to their high electrical/thermal conductivity or flame retardancy. However, to fully employ the benefits of polymer/EG nanocomposites one must consider the high degree of porosity of EG. The high degree of porosity of EG can deteriorate the composites’ mechanical properties if the polymer chains cannot diffuse completely into the EG pores. In this article, an insight is given into the diffusion of unsaturated isophthalic polyester (UP) resin, consisting of a combination of maleic anhydride and isophthalic anhydride in the resin backbone, with two viscosities, into the pores of the EG particles of various degrees of porosity. The diffusion experiments were carried out on compressed EG tablets with the same density but different porosity due to the different porosity of the EG particles. The results showed that the diffusion rate of the UP resin with higher viscosity slightly decreased when the EG porosity decreased but, in the opposite way, it strongly increased for the low viscosity UP resin. The EG nanocomposites samples were molded at varying pressures. The micrographs of the fractured surfaces of the EG nanocomposites showed that the EG pores were not filled with resin, thus the EG nanocomposites had residual pores. It was found that composites containing EGs with higher expansion ratio and larger particles and pores showed larger residual pores. Furthermore, the composites prepared with the more viscous UP resin showed more residual pores. By applying a pressure of 10?bar instead of 1?bar, a reduction of 7–20% in the residual pores of the nanocomposites was observed which led to improved mechanical properties by up to 20% in flexural strength for the EG with the highest expansion ratio.     

Monte Carlo simulations of magnetic particle suspensions with a simple assessment method for the particle overlap between magnetic spheroids

We have developed a simple assessment method for the overlap between spheroidal particles, which neither requires the complex manipulation of vectors and matrices that is indispensable in the ordinary methods, nor is based on a model potential. Moreover, we have developed an evaluation method for the interaction energy arising from the overlap of the steric layer coating spheroidal particles. This is based on a sphere-connected particle model, but some modifications are introduced in order to express an appropriate repulsive interaction energy at the deepest overlapping position. We have investigated the phase change in a magnetic spheroidal particle suspension for a two-dimensional system by means of Monte Carlo simulations. In the case of no external magnetic field, if the magnetic particle-particle interaction is sufficiently strong to favour cluster formation, long raft-like clusters tend to be formed in a dilute situation. With decreasing values of area fraction, a chain-like structure in a dense situation transforms into a raft-like structure within a narrow range of the particle area fraction. Similarly, the raft-like clusters are preferred in a weak applied magnetic field, but an increase in the field strength induces a phase change from a raft-like into a chain-like structure.

Highlights of the present paper:

1. A simple assessment method has been proposed for the overlap between two spheroidal particles.

2. The particle overlap assessment is free from a complex mathematical manipulation regarding vectors and matrices.

3. A modified sphere-connected model has been proposed in order to more accurately evaluate a repulsive interaction due to the overlap of the steric layers coating spheroidal particles.

4. 2D Monte Carlo simulations have been performed to elucidate the phenomenon of a phase change by magnetic spheroidal particles on a material plane surface.

5. A phase change between a raft-like and a chain-like aggregate structure is able to be controlled by the area fraction of particles and an external magnetic field.

     

利用电磁法研究HMX与TATB混合钝感炸药的冲击起爆特性

为研究含有少量奥克托金(HMX)且以三氨基三硝基苯(TATB)为基的高能钝感炸药PBX-3的冲击起爆反应增长规律,采用火炮驱动蓝宝石飞片的方法和铝基组合式电磁粒子速度计技术进行了一维平面冲击实验。通过实验测量撞击表面及内部不同深度处的冲击波后粒子速度,得到PBX-3炸药的Hugoniot关系。根据冲击波示踪器所测数据绘制了炸药到爆轰的时间-距离(x-t)图,获得了反映炸药冲击起爆性能的Pop关系。将入射压力为12.964 GPa时达到爆轰的6条速度曲线修整成相同零点,通过读取6条曲线的分离点即反应区末端的C-J点,计算出化学反应区时间和宽度。     

运用三维离散元技术模拟落锤撞击下奥克托今颗粒的点火燃烧过程

炸药颗粒的点火燃烧过程一直是人们关注的热点问题。近年来，三维离散元技术在中尺度观测颗粒材料的动力学过程中拥有显著优势。炸药燃烧属于颗粒材料的反应动力学，运用三维离散元技术(DM3)可以有效地观测炸药燃烧传播的过程。以奥克托今(HMX)颗粒为例，本文成功模拟并观测到了HMX颗粒的燃烧反应程度，确定了颗粒开始燃烧反应的时间，以及燃烧反应传播的时间。同时，结合落锤冲击颗粒的三维图像以及其表观压强和放热功率，得到了HMX颗粒燃烧反应、燃烧传播的整个反应动力学过程，包括颗粒在冲击加载下碎化塑性变形的过程，颗粒燃烧反应放热的过程，落锤回弹颗粒喷射的过程等。同时，进一步说明了尖顶颗粒更利于颗粒点火，平顶颗粒有抑制颗粒点火的能力。